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We examine the global structure of scalar field critical collapse spacetimes using a characteristic 
double-null code. It can integrate past the horizon without any coordinate problems, due to the 
careful choice of constraint equations used in the evolution. The limiting sequence of sub- and 
supercritical spacetimes presents an apparent paradox in the expected Penrose diagrams, which we 
address in this paper. We argue that the limiting spacetime converges pointwise to a unique limit 
for all r > 0, but not uniformly. The r = line is difi'erent in the two limits. We interpret that the 
two different Penrose diagrams diff'er by a discontinuous gauge transformation. We conclude that 
the limiting spacetime possesses a singular event, with a future removable naked singularity. 
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I. INTRODUCTION 

The critical collapse of scalar fields gives rise to a new 
class of thought experiments in general relativity P, 01 ■ 
It has been suggested that a weak singularity may be vis- 
ible to a distant observer during the collapse of a scalar 
field Si. With the presence of such potentially prob- 
lematic phenomena, one can ask to what degree cosmic 
censorship is violated. The spirit of the cosmic censorship 
conjecture is that the evolution of spacetimes as seen by 
distant observers in asymptotically flat regions is deter- 
mined uniquely by the Einstein equations. If visible infi- 
nite curvature arises in the evolution from regular initial 
data, one might have cause for concern about the classical 
completeness of Einstein's theory of gravity. The critical 
collapse is a fine-tuned limit, and "strong cosmic cen- 
sorship" has been formulated to exclude such rare cases 
with zero measure in the space of initial conditions. It is 
nevertheless instructive to understand the nature of this 
singularity, and how such a limit can be taken. 

Recently, Martin-Garcia and Gundlach have numer- 
ically constructed a self-consistent discretely self-similar 
spacetime which they argued to be related to the evo- 
lutionary critical collapse solutions. They proposed the 
existence of a future Cauchy horizon emanating from the 
critical collapse event, on which new data must be spec- 
ified. The authors proceeded to find a unique way of 
specifying this data which can result in a regular future 
spacetime. 

In this paper, we consider the problem from a different 
perspective. Instead of searching for the critical solution 
of Einstein equations by first taking the limit of discrete 
self-similarity, we study the constructive sequence of non- 
critical global spacetimes. We then search for a unique 
limit as we approach criticality. Posed in such a way, the 
existence of a Cauchy horizon would be very unexpected: 
each spacetime in the limit sequence has no Cauchy hori- 
zon, so why would it form in the limit? To study the 
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problem, we developed a characteristic code which can 
track a scalar field collapse interior to the horizon. In Sec- 
tionm we outline a conceptual paradox in the search for 
a global critical spacetime structure. In Section ITTll we 
derive the scalar field equations of motion with spherical 
symmetry. We proceed to solve these equations numeri- 
cally in Section llVI The numerical results are presented 
in Section where we explain our proposed solution to 
the apparent paradox. 



II. THE APPARENT PARADOX 

Consider a collapsing scalar field with amplitude char- 
acterized by a parameter p. For small amplitudes p < p<,, 
the field collapses and re-expands. For large amplitudes 
p > p^,, a black hole forms. An interesting question is to 
examine the behavior when p — p*. One can consider two 
limits, one from below and one from above. If we con- 
sider a sequence of subcritical spacetimes with collapsing 
scalar fields as they approach p^ from below, we might 
expect the global structure of the resulting spacetime to 
be Minkowski, i.e. a triangle shown in the left panel of 
Fig. ^ If we consider the limit of a sequence of supercrit- 
ical spacetimes from above, we expect each stage to have 
a global Schwarzschild structure. As the parameter de- 
creases, the mass of the resulting black hole approaches 
zero, and one expects the limiting spacetime to resemble 
a Schwarzschild spacetime with zero mass, shown in the 
right panel of Fig.^ In this limit, r = coincides with 
the horizon and becomes null (and infinitely redshifted). 
The global structure looks quadrangular, quite unlike the 
argued spacetime in the other limit. 

This suggests several possible interpretations. Perhaps 
the two limits are different, and the limiting spacetime 
depends on the direction from which the limit was taken. 
Or one of the limits is only an incomplete description of 
spacetime, and might be extensible to the same global 
structure. Or perhaps the limit is not convergent from 
either direction, and oscillates in such a way that more 
data must be specified on a spontaneously formed Cauchy 
horizon . Our study suggests a slightly different phys- 
ical interpretation of the global structure. 
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FIG. 1: Two possible global structures of a critical collapse 
spacetime. The left is the diagram expected taking a limit of 
subcritical collapses, while the right diagram is a zero mass 
black hole resulting from the limit sequence of supercritical 
coUapse. 



Critical spacetimes are only known as numerical solu- 
tions, which makes questions about global structure hard 
to answer. It is most easily studied in characteristic co- 
ordinates which follow light ray propagation '6]. In the 
subsequent section we will describe our formulation and 
implementation of the numerical procedures. 



III. SPHERICALLY SYMMETRIC SCALAR 
FIELD COLLAPSE 

The spherically symmetric (n -I- 2)-dimensional space- 
time metric can be written as 



(1) 



where dfl^ is the metric of a unit 7i-dimensional sphere 
with curvature K — 1, and the two-manifold metric 



^dx is conformally flat. The dynamics of the 



scalar field collapse are described by the reduced action 
S (X /d2a;e"''|-2(V0)2-2n(V(7-VM) (2) 



where integration and differential operators are with re- 
spect to the fiat two- metric dx^. Variation of the above 
action with respect to the fields </>, ^ and cr gives equa- 
tions of motion 

□0 + n{Vii ■ Vcj}) = 0, (3a) 
Ufi + n{Viif - {n - l)i^e-2('^+M) = q, (3b) 
□(T-§{nAi+(VAi)2}-(V0)2 = O, (3c) 

while the two constraint equations are recovered by vari- 
ation with respect to the (flat) metric 

Traceless {/i;ah + M.aM,b + 2Ai(,aO-,&) + f 0,a<^,&} = 0. (4) 



For the particular case of spherically symmetric scalar 
field collapse in four dimensions (n — 2), the equations 
of motion Q can be simplified by introducing auxiliary 
field variables r — and — r(j): 



□ (r2) = 2e-2'^, 



(5a) 
(5b) 
(5c) 



This form of the dynamical equations is better suited for 
numerical evolution. 



IV. CHARACTERISTIC CODE 

We discretize and evolve the collapsing scalar field 
spacetime in double- null coordinates dx^ — ~2dudv, 
where the radial characteristics of the wave equations 
are made explicit: the outgoing characteristic propagates 
at constant u in the direction of increasing v, while the 
incoming characteristic propagates at constant v in the 
direction of increasing u. This approach has a number of 
advantages over some of the more traditional spacetime 
slicings. 

The characteristic code only propagates information 
along characteristics at a numerical speed equal to the 
true characteristic speed. The numerical domain of de- 
pendence is no larger than the physical domain of depen- 
dence, and the code still maintains full (second order) 
accuracy. Horizons are not particularly special as far as 
ingoing null characteristics are concerned. This allows us 
to follow collapse all the way to the singularity. Even as 
an outgoing characteristic hits a singularity and the float- 
ing point numbers denormalize, this does not affect any 
of the other characteristics, which can be integrated to 
fill the whole maximally extended spacetime determined 
by the initial data. To illustrate this point, in Fig. [3 
we present Penrose and Kruskal diagrams of a spacetime 
with a large black hole formed in the collapse of the scalar 
field wavepacket. 

The two-dimensional metric c?7^ has a residual gauge 
freedom under redefinition of the null coordinates u i— > 
U{u), V 1-^ V{v). These two free functions (of a single 
variable) are used to define the coordinate v on the initial 
slice and to map the central point r = to a straight line 
in the (it, v) plane: 



v\u=o^V2r, u\r=o = v. 



(6) 



The second gauge condition is particularly convenient, 
since it places the central point at a known location on 
the grid when discretizing. 

The initial conditions are specified on a surface of con- 
stant u by giving a scalar field profile (j>{v). Together 
with the gauge choice lO, this determines the rest of the 
variables. In particular, cr(w) is obtained by integrating 
the outgoing (vv) constraint equation Q 



a\u=o i (l^^vdv. 



(7) 
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FIG. 2: Penrose (left) and Kruskal (right) style diagrams of the spacetime with a large black hole produced in the collapse of 
a sine-squared scalar field wavepacket (I14II . The grid lines on the Penrose diagram show observers at constant radius and the 
proper time according to their clocks (note that the plot is rotated 45° compared to Fig. Q the diagonal line u = v is the origin 
r = 0). The Kruskal diagram shows the trajectories of outgoing null rays as emitted from the center of spacetime at subsequent 
moments of time. The code sees both real and apparent horizons, as well as the formation of the spacelike singularity (small 
region around which is excised). Note that the apparent and event horizons are not the same, as the spacetime is not static. 



The integration is implemented as a fourth order Runge- 
Kutta algorithm with a fixed step. Although strictly 
speaking it is not necessary, as the evolution code is sec- 
ond order, it is no more complicated than a second order 
integrator would be. 

Although the constraint equations follow from the dy- 
namical equations, their use might be required for stabil- 
ity 0- Using the outgoing constraint equation for evo- 
lution (rather than just for initial conditions) is not a 
good idea, however. It becomes degenerate on the ap- 
parent horizon, where the outgoing light rays become 
(marginally) trapped ~ 0. The incoming constraint 
equation is perfectly fine, though, and could be followed 
all the way to the center of the spacetime, whether it is 
singular or not. By introducing an auxiliary field variable 
T — CT-I- i In (^—^/2r,u), the incoming constraint equation 
Q can be written in a form that is very simple to inte- 
grate 



tion of the evolution equations. The covariant differential 
operators in the evolution equations (jSJ are with respect 
to a flat metric, so in null coordinates they are written 
simply as 

Dx = -2dudyx, (Vx)^ = -2{dux){dyx), (10) 

where x denotes any one of the three dynamical variables 
in equations (O. We discretize by finite differencing on 
a rectilinear (m, v) grid with equal spacing du = dv — A, 
which to second order accuracy gives 



2 

(□a;)x = -^[^^ 



x+- - x-+\ , 



(11) 



(Vx) 



2A2 



[{x++ - X — f - {x-+ - X+-)'^ 



where the differential operators are evaluated at the cen- 
ter of a grid cell (see Fig. OJ. The code takes a step by 



2(lnr),„ 



(8) 



To integrate the incoming constraint, one would need to 
know the values of r on the initial slice u = const. These 
can be found by integration of the equation 



1 ,2 1 r 



-2rl 



(9) 



which is obtained by combining the constraints with the 
evolution equation (|S^) to solve for on the initial slice. 

Having discussed our gauge choice, initial conditions 
and constraint equations, we now come to the discretiza- 



X 



FIG. 3: Discretization on a (m, v) grid. 
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using discretized evolution equations to find the values of 
the fields at the (++) grid point (shaded node in Fig.|2J). 
The only non-trivial operation involved in this is finding 
r_|_+ accurate to second order, which is done by solving 
equation jsji) discretized in the following fashion: 



2 2 2 

r — r I — r_ 



A , 
^exp(-r. 



+ ^ 



(12) 



The rest is then straightforward, as the right hand sides 
of equations (03, c) become known. 

As the code advances to the next slice of constant u, 
the very first point on the grid (u = v) is the center of 
the spacetime (r = 0) and has to be treated specially. 
The asymptotic form of the evolution equations at r = 

Vr-V0 = O, cr = -iln(Vr)2 (13) 
is used to calculate the values of the fields in the center. 

V. RESULTS 

We tested the code on the collapse of a scalar field with 
various initial conditions. In particular, we used pulse 



_ J psin 4:TTv, -J < V < 



otherwise 



and kink 



{v) = 



ptanh [5 tan7r(|i; — i)] , v< 
2p, V > 



(14) 



(15) 



field profiles in our simulations. Both profiles are fairly 
smooth functions with field energy density having com- 
pact support on the initial slice. This avoids the interfer- 
ence of long-range tails in the initial data on the late-time 
evolution. 

All runs shown in this paper used 65536 uniformly 
spaced grid points. The critical scaling of black hole mass 
is shown in Fig.0| It agrees well with the literature 
Our achievable dynamic range shows scaling over three 
orders of magnitude in black hole mass without the use 
of adaptive mesh refinement. This is sufficient for our 
study. Further improvements in dynamic range can be 
sought after using adaptive mesh refinement techniques 
or adopting an initial gauge which places ingoing null 
rays on a grid more densely near the collapse point ^^j. 
As the grid resolution is increased, the double arithmetic 
precision and round-off errors become the main obstacles 
in the quest for higher dynamical range. 

Fig.Elshows the Penrose diagram of a subcritical space- 
time just below the threshold of black hole formation. 
The lines of constant radius r are drawn as solid curves, 
and the proper time along such lines are drawn as dashed 
curves. Most of the scalar field mass is shed before 



-11 -10 -9 



-4 -3 -2 -1 



ioaio(P - P-) 

FIG. 4: Black hole mass scaling for single kink wave, given 
by equation 11511 . The upper panel shows a power law fit to 
the numerical data with critical exponent 7 ~ 0.373 which 
holds over three orders of magnitude in black hole masses, 
corresponding to nine orders of magnitude in the tuning pa- 
rameter. The lower panel shows the residual after the power 
law fit, demonstrating periodic fine structure due to discrete 
self-similarity of the collapse. Its period lu « 4.63 corresponds 
to the value of the echoing parameter A « 3.45. The lower 
panel inset shows the residual mapped to a single period (out- 
lier points removed) , illustrating that the numerical data is in 
good agreement with a periodic modulation of the black hole 
mass. 



u ~ 0.6, but the field oscillating on ever smaller scales 
creates near-singular curvature in the center of space- 
time. 

When we examine the Penrose diagram of a supercrit- 
ical collapse as shown in Fig. IHl we see the formation of a 
horizon near u « 0.55. The code evolves well inside the 
event horizon, and we can identify the apparent horizon 
in the interior of the black hole. The lines of constant r 
become spacelike beyond the apparent horizon. All the 
larger radii outside the black hole pile up at the horizon 
in this diagram. 

We can pose the question if the spacetime described 
by Fig. could possibly be the same as Fig. [S] Our 
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FIG. 5: Penrose diagram for subcritical collapse. 



r = const : 

I = const — 7 - - 
apparent tiorizon - 

event horizon . 



0.8 - 



0.6 - 




0.2 0.4 0.6 0.8 1 



FIG. 6: Penrose diagram for supercritical collapse with 
Mbh/Mo « 0.046. 
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FIG. 7: Re-gauged supercritical collapse of Fig. |^ The re- 
gauging was done at the postprocessing stage from the same 
simulation run. One sees the clear similarity with the sub- 
critical collapse shown in Fig. |3 



original gauge choice from equation © fixes the r = 
line to coincide with u — v. It does not have to be 
this way; the null coordinates leave the possibility for 
a global gauge change u i— > U{u). One can try to "un- 
pile" the lines of constant r observers for the supercritical 
collapse spacetime in Fig. El by defining a gauge change 
U{u) such that the second r curve coincides with that in 
a subcritical collapse spacetime of Fig. [S] The result is 
shown in Fig. |7| 

This stretching is rather sudden and throws the r = 
line into an almost null direction, but all the other fi- 
nite r's appear to coincide. This is non-trivial, since the 
residual gauge freedom is a single function of one vari- 
able, which we used to fix one value r. If the spacetimes 
are equivalent, the other r and t should fall into place, 
and similarly if the spacetimes are different they should 
diverge. We indeed see that the proper time and outer 
radii move into place, as would be expected in a conver- 
gent spacetime. It is also clear that the r = line is not 
convergent, and the limiting spacetime experiences infi- 
nite curvature at that one point. To address the future 
of this singular event, we can look at the spacetime of 
the limiting sequence of super- and subcritical collapses. 

We now discuss the global structure of the critical 
spacetime. Our first question is the nature of the r = 
line in the course of the singular collapse event. An ob- 
server moving at some finite radius r sees the collapse of 
a scalar field to a point, and a rebound. This is true for 
both the supercritical and subcritical collapse, since even 
when a black hole forms, the majority of the initial scalar 
field energy escapes, and the black hole only contains an 
ever smaller fraction of the initial mass as the parameter 
is tuned to criticality. Long after the field rebounds, the 
observer can measure the gravitational redshift to adja- 
cent interior radii, and concludes that there is no red- 
shift for most radii. In the supercritical scenario, there 
are significant redshifts at scales r < Mbh- But for any 
fixed radius observer, as one takes the limit of p — > p*, 
the sphere of influence of the ever diminishing black hole 
mass shrinks to zero, and the spacetime converges point- 
wise to Minkowski at r > 0. The convergence is not uni- 
form, since for any given redshift difference e between the 
observer and a fiducial interior radius, one can construct 
a radius r M/e inside which the redshift is larger than 
e. In the two limits p p*^_ and p ~* p*_, the spacetimes 
converge pointwise everywhere except for the line r = 0. 

A simple analogy is a singular weak field star. Con- 
sider a sequence of spacetimes with a single star of mass 
and radius = A//*/e. The Newtonian potential out- 
side the star is V{r) — —M^jr for r > M^jt. To simplify 
the argument, we make the potential continuous and con- 
stant interior to its radius V(r) — e. For small values of 
e, the spacetime is in the weak field regime everywhere 
to order e. If we take a sequence of such space times 
with fixed e and decreasing M^, the maximal curvature 
R cx €^ /Ml increases without bound. There appears to 
be an illusory naked singularity in that limit. We can ask 
what the limiting spacetime looks like as M, — > 0. The 
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obvious answer would be empty Minkowski space. The 
convergence to this space is pointwise, but not uniform. 
In the limit, Vo{r) = limM^^o M^/r. Vo{r) converges to 
pointwise, but not uniformly, so in the limit, Vo(r) = 
everywhere except for r — 0, and it is undefined at that 
point. This is a removable singularity since we can define 
Vo(0) = limr.^0 = 0, which is a unique and physi- 

cally acceptable solution. We suggest that the future of 
the singular collapse is likewise regular everywhere, with 
a removable singularity at r = in the future of the 
collapse event. 

The apparent paradox in the Penrose diagram seen in 
Fig.CJarises from the gauge choice at r = 0. This one line 
is poorly defined in the limit, since it becomes singular. 
The left panel describes the physical spacetime in the 
critical collapse limit, and the right one is related by a 
singular gauge change. 

Our analysis takes a very different approach from 
Martin-Garcia and Gundlach 5] . These authors solve an 
elliptic equation satisfying an ansatz of discretely self- 
similar critical collapse. We took a limit of a series of 
hyperbolic initial value problems. Those authors found a 
possibility of specifying new Cauchy data in the elliptic 
solution, which is never an option for our evolutionary 
approach. The qualitative solution for their unique reg- 
ular extension looks similar to our critical spacetime. 

VI. CONCLUSIONS 

We have implemented a characteristic code in double- 
null coordinates and used it to study the global spacetime 



structure of a critical collapse of a scalar field. We repro- 
duce the standard critical behavior and universal scaling. 
We compare the limiting spacetime from the subcritical 
and the supercritical collapse limits. These two limits 
appear qualitatively different. Based on the numerical 
simulations, we conjecture that the two limits converge 
pointwise to the same spacetime, but not uniformly. In 
particular, the r = line is not convergent, but all other 
points appear to converge pointwise. The two apparently 
different solutions then only differ by a gauge change. 
The apparent naked singularity in the upper limit ob- 
tained by the sequence of spacetimes with ever decreas- 
ing black hole mass becomes a removable singularity in 
the limit. 

We conclude that the nature of the limiting critical 
collapse space time requires a careful definition of the 
order in which limits are taken, since the convergence is 
not uniform. For any collapse parameter p within the 
critical value \p — p*\ = S > 0, one can find sufficiently 
small radii within which the supercritical and subcritical 
solutions differ. Conversely, for any fixed r > 0, one can 
find a value 6{r) where the super and sub critical space 
times agree to some tolerance e at radius r. There is a 
unique limit, with no new Cauchy data that emits from 
the singular event. 
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